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Green’s function and anti-holomorphic 
dynamics on a torus 

Walter Bergweiler and Alexandre Eremenko* 


Abstract 

We give a new, simple proof of the fact recently discovered by C.- 
S. Lin and C.-L. Wang that the Green function of a torus has either 
three or five critical points, depending on the modulus of the torus. 
The proof uses anti-holomorphic dynamics. As a byproduct we find 
a one-parametric family of anti-holomorphic dynamical systems for 
which the parameter space consists only of hyperbolic components 
and analytic curves separating them. 


1 Introduction 


Green’s function on a torus T is defined as a solution of the equation 


AG 


-5 + 


1 


in’ 


normalized so that 



0 . 


Here S is the delta-function, and |T| is the area of T with respect to a flat 
metric. 

We write the torus T as T = C/A with a lattice 


A = {mui + nu2: m,n E Z}, 
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where r = UJ 2 /UJ 1 satisfies Imr > 0. Recently C.-S. Lin and C.-L. Wang [ 21 ] 
discovered that Green’s function has either three or five critical points, de¬ 
pending on r. It is surprising that this simple fact was not known until 2010. 
In [H [To] they study the corresponding partition of the r-half-plane. Their 
proofs are long and indirect, using advanced non-linear PDE theory. Our 
paper is motivated by the desire to give a simple proof of their result that 
Green’s function has either three or five critical points and to give a criterion 
for r distinguishing which case occurs. 

We have (see m) 


G{z) 


1 


log 101 (z) I + 


(Im^;)^ 
2 Imr 


+ C'(t), 


where 61 is the first theta-function. Here and in the following we use the 
notation of elliptic functions as given in mm- We note that the notation 
in [H [32] [33] is different, see the remark following the theorem below. 
Gritical points of G are solutions of the equation 


C{z) az + bz = 0 , 


( 1 ) 


where the constants a and b are uniquely defined by the condition that the 
left hand side is A-periodic. With ({z + Uj) = ({z) + rjj for j = 1, 2 we thus 
have 


rji + aui + bui = 0 and 772 + cb ^2 + buj 2 = 0. 
With the Legendre relation r]iUJ 2 — r] 2 <^i = 27rz we obtain 


b = 


TT 


IcuiP Imr 


and a = 


buji rji 


TT 


Ui 


Vi 


ixJi uf Im r (Ui 


So the problem is to determine the number of solutions of ([^ where 
given by (|^. 


are 


Theorem. The equation 0 has three solutions in T if CjOjl -1- rjiOJi 


( 2 ) 

a and b 

= 0 or 


Im 


2Tri 

ejuj + r]iui 



> 0 


( 3 ) 


for some j G {1, 2, 3} and it has five solutions otherwise. 
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Here, as usual, ei = p{uJi/2), 62 = p{{oJi +oj 2 )/ 2 ) and 63 = p{pj 2 / 2 ). An 
elementary computation shows that the condition in the theorem is equivalent 
to 

min 

i<i<3 

We note that CjCnf and r^icui depend only on r = We may restrict 

to the case that Wi = 1 so that t = 102 - Then ([^ simplihes to 

27ri \ 

-r > 0. 

ej + ?7i J 

As mentioned, a different notation for elliptic functions is used in [21 [321 
l33] . There the periods are denoted by 2uj and the definition of rji also differs 
by a factor 2. Thus in that terminology (|^ takes the form 

ni 

2 {ejul + Tjiuii) 

Figure [^shows (in gray) the regions in the r-plane where Green’s function 
has 5 critical points; that is, the set of r-values where (|^ fails for all j. The 
range shown is | Rer| < 1 and 0.15 < Imr < 2.15. 

The standard fundamental domain consisting of those r which satisfy the 
inequalities Imr > 0, —1/2 < Rer < 1/2 and |r| > 1, with |r| = 1 only if 
Re r > 0, is in the upper middle of the picture. Its images under the modular 
group are also shown. 

Our proof is based on Fatou’s theorem from complex dynamics. Orig¬ 
inally Fatou’s theorem was proved to estimate from above the number of 
attracting cycles of a rational function. Then it was extended to more gen¬ 
eral classes of functions. The most surprising fact is that Fatou’s theorem 
can be used sometimes to estimate the number of solutions of equations in 
settings where dynamics is not present. This was first noticed in [12]; the 
contents of this unpublished preprint is reproduced in [T31IH|- The paper [U] 
shows that Fatou’s theorem can be used to prove under some circumstances 
the existence of critical points of a meromorphic function. In the papers 
[19] Fatou’s theorem was used to obtain upper estimates of the 
numbers of solutions of equations of the form 

= r(z), 



r > 0. 



ejojf + rjiOJi 


Imr — 1 


TT 


< 1 
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Figure 1; The regions given by (|^ in the r-plane. 


with a meromorphic function r; this permitted to prove a conjecture in as¬ 
tronomy | 2 ni El]- In the recent work [23|, a topological classification of 
quadrature domains is obtained with a method based on Fatou’s theorem. 

As mentioned, Fatou [T^ §30] stated his result originally only for rational 
functions. As pointed out for example in |5l Theorem 7], the proof extends 
to functions meromorphic in C. 

To state the version of Fatou’s theorem that we need, we recall some defi¬ 
nitions. Let S and T be Riemann surfaces and let /: S' —)■ T be holomorphic. 
A point c G S' is called a critical if f'{c) = 0; this condition does not depend 
on the local coordinates. 

We say that a curve 7 : [0,1) —)■ S escapes if for every compact subset 
R in S' there exists to ^ (0,1) such that 'jit) ^ K for t E [to,l)- A curve 
7 : [0,1) —)■ S' is called an asymptotic curve of / if 7 escapes and the limit 
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limt_j.i /( 7 (t)) exists and is contained in T. This limit is called an asymptotic 
value of /. The following result is due to Hurwitz na. 

Hurwitz’s Theorem. If a holomorphic map between two Riemann surfaces 
has no critical points and no asymptotic curves, then it is a covering. 

This result can be deduced from the fact that covering maps are charac¬ 
terized by the path-lifting property: for every curve 7 : [0,1] —?• T and every 
Zq G /~^( 7 ( 0 )) there exists a unique curve 7 *: [0,1] —)• S' such that 7 = / 07 *; 
see, for example, |21 Theorem 9.1], where the term “complete covering” is 
also used. Note also that critical points and asymptotic curves correspond 
to singularities of the inverse function; see [2HI §XI.l]. 

A Riemann surface S is called hyperbolic if its universal covering is the 
unit disk [2] . A hyperbolic Riemann surface is equipped with the hyperbolic 
metric. We denote by Xs{z)\dz\ the length element of the hyperbolic met¬ 
ric in S. The invariant form of the Schwarz lemma (see, for example, [STl 
Proposition 1 . 2.8 (a)], m Theorem 2.11] or |H Theorem 10.5]) says that a 
holomorphic map /: S' —?• T between hyperbolic Riemann surfaces S' and T 
satishes XT{f{z))\f{z)\ < Xs{z), with strict inequality for all z, unless / is 
a covering, in which case the equality holds for all z. 

For a holomorphic map /: S' —)■ S' a point Zq ^ “S' is called fixed if /(^o) = 
zq, and for such a point f'{zo) is called the multiplier. It is easy to verify 
that the multiplier does not depend on the local coordinate. A hxed point is 
called attracting, neutral or repelling depending on whether the modulus of 
its multiplier is less than, equal to or greater than 1 , respectively. 

Fatou’s Theorem. Let S be a hyperbolic Riemann surface and let f: S ^ S 
be a holomorphic map with an attracting fixed point zq E S. Then f has 
a critical point or an asymptotic curve in S'. Moreover, /"■ —)• zq locally 
uniformly in S'. 

Proof. If / has no critical points and no asymptotic curves, then / is a cov¬ 
ering. The Schwarz lemma yields As( 2 ;o) = A 5 (/(^o))|/'(^o)| = Xs{zo)\f'{zQ)\ 
and thus \ f\zo)\ = 1, contradicting the assumption that \f'{zo)\ < 1. The 
second statement also follows from the Schwarz lemma. □ 

In applications to holomorphic dynamics, S is the immediate attraction 
basin of an attracting hxed point Zq] that is, the component of the Fatou 
set which contains Zq. It then follows that the immediate attracting basin 
contains a critical point or an asymptotic curve. 
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2 Proof of the Theorem 

To prove the theorem we rewrite ([^ as a hxed point equation 

z =-^(^C{z)+ =: g{z), (4) 

so g is an antz-meromorphic function in the plane. As the left hand side 
of Q is A-periodic, we conclude from (|^ that 

g{z + u) = g{z) + u, cu G A. (5) 

The function g does not map the plane into itself because it has poles, 
so the equation (|^ does not permit to dehne a map of the torus T into 
itself. To remedy this, we consider the set Pq of poles of g. For n > 1 we 
dehne inductively P„ = g~^{Pn-i). Then all iterates of g are dehned on 
the set C\Poo, where Poo = Pn- Let J be the closure of Poo. Then 
the iterates of g form a normal family in C\J, and this set is completely 
invariant under g. Thus - with an obvious extension of these concepts from 
holomorphic functions to anti-holomorphic ones - we call J the Julia set 
and its complement F = C\J the Patou set of g. The Fatou set is thus 
the maximal open subset of the plane such that g{F) C P. Evidently, P is 
A-invariant, so the map g: F ^ F descends to a map which is dehned on an 
open subset of the torus T and maps this open subset to itself. 

Let vr: C —?• C/A = P be the projection map, P = vr(P), and g: F ^ F 
the induced map which satishes 'g o n = nog. The Riemann surface P is a 
subset of the torus, and it is hyperbolic because its complement is inhnite. 

We apply the terminology attracting, neutral and repelling also to hxed 
points Zq of anti-holomorphic maps /, considering df{zo) as the multiplier. 
Notice that 

9giz) = -^iCiz) + a) = ^{p{z)-a). (6) 

To obtain holomorphic dynamics instead of the anti-holomorphic one, we 
consider the second iterates h = g^ and h = ^. Then we have 

h O 71 = 71 o h. ( 7 ) 

Images of hxed points of g and h in P under tt are hxed points of g and h, 
respectively. By the chain rule, we have 

h'= {(dg)og)-dg = {(dg)og)-dg. (8) 
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For a fixed point Zq of g we thus obtain h'{zo) = \dg{zQ)\‘^. Even though we 
will not need this fact, we observe that the multiplier of h at a hxed point of 
g is always a non-negative real number. 

In order to apply Fatou’s theorem to h: F ^ F, we have to consider the 
critical points and asymptotic curves of h. 

Lemma 1. The map h: F ^ F has no asymptotic curves. 

Proof. We prove this by contradiction. Let 7 be an asymptotic curve of h 
in F. Let 7 be some lifting of 7 in F] that is tt o 7 = 7 . As 7 escapes 
from F, the curve 7 escapes from F. By assumption, we have h{j{t\) —>■ p 
as f —)■ 1, for some p G F, so h{ 7 r{'y{t))) —)■ p as f —)■ 1. In view of ([7) this 
yields 7 i{h{'y{t))) —)■ p as t —>■ 1. Because vr is a covering we conclude that 

h{'y{t)) —)■ p for some p E F, (9) 

with 7 r(p) = p. Thus 7 is an asymptotic curve of h in F. 

Suppose that 7 has a limit point g G C as t —>■ 1; that is, there exists 
a sequence (tj) tending to 1 such that 7 (tj) —)■ q. Then q G dF C J. If 
q is neither a pole of g nor a preimage of such a pole under p, then h is 
holomorphic at g, with h{q) G J by the complete invariance of the Julia set. 
On the other hand, we have h{q) = limj_,.oo h{'j{tj)) = p E F, which is a 
contradiction. If q is the preimage of a pole of g, then g is a pole of h. Thus 
h(g) = 00 , contradicting again h(g) = p. 

This shows that the finite limit points of 7 are poles of g. As the poles 
form a discrete subset of C we actually see that if 7 has a hnite limit point 
g, then 'jit) —)■ g as t —)• 1 and thus a(t) := g^'jit)) —>■ cx) as t —)■ 1. On 
the other hand, g{a{t)) = h{'j{t)) —)■ p as t —)■ 1. Let F be a A-invariant set 
consisting of disjoint disks around the poles of g. It follows from (|^ that 

g{z) -E 00 as z —>■ 00 , z G C\F. 

This is a contradiction to a(t) —)■ 00 and g{a{t)) —)■ p as t —)■ 1. We have thus 
shown that 7 has no hnite limit points, meaning that 7 (t) —)■ 00 as t —?• 1 . 

In order to show that this is impossible we note again that the singularities 
of h are the poles of g and their preimages under g. Since these preimages 
accumulate only at the poles of g (and at cxd) there exists a A-invariant set 
F' C C consisting of disjoint disks which contains all singularities of h. Since 
h satishes 

hi^z J- (u) = h(yZ^ J- (u, uj G A, 
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by 0 , it follows that 

h{z) —)■ oo as z —)■ cx), z G C\E'. 

As before this is incompatible with 7 (t) —)■ oo and —)■ p as t —)• 1. 

This completes the proof of the lemma. □ 

We consider the eqnivalence relation on F defined by z ~ z' if h'^{z) = 
h^i^z') for some non-negative integers m and n. The eqnivalence classes are 
called the grand orbits. We note that if the seqnence {h'^{z)) converges for 
some z E F, then for all z' in the grand orbit of ^ the sequence {h'^{z')) 
converges to the same limit. We call this the limit of the grand orbit. 

Lemma 2. The set of critical points of h belongs to at most 4 grand orbits 
under h. 

Proof. By (|^, the equation dg{z) = 0 is equivalent to p(z) = a, and has 
two solutions modulo A which dehne points c and —c on T. Here and in 
other similar places we denote by —c the point which corresponds to c by the 
conformal involution of the torus. If any of the points c and —c is in F, then 
it is a zero of h'. By ([^ the other zeros of h' are p“^(±c) = ±p“^(c). Even 
though this is an inhnite set on T, the critical points of h' are thus contained 
in at most 4 grand orbits represented by c, —c, p(c) and —p(c). □ 

Lemma 3. Let Zq be a fixed point ofg. If h'^{z) -E- Zq as n ^ oo for some z, 
then hP(g{z)) -E zq. 

Proof. We have h"' o If ='g o h'^ for all n G N. Since h"‘{z) -E Zq this yields 
= g(jp{z)) -E g{zQ) = zq. □ 

Applying Eaton’s theorem to the map h: F ^ F, we deduce that h and 
thus g has at most two attracting hxed points. If F is disconnected, we apply 
Eaton’s theorem separately to each component that contains an attracting 
hxed point. As the vr-image of an attracting hxed point of g is an attracting 
hxed point of 'g we obtain 

Lemma 4. The function g has at most two attracting fixed points, modulo A. 

The map 0: T —>■ C, 2 ; i—)■ z—g{z), is well dehned by ([^. Let = l — \dg\^ 
be the Jacobian determinant of 0. Then = {z E T\ J(j,{z) > 0} is the 


set where 0 preserves the orientation and D~ = {z E T: J^{z) < 0} is the 
set where (j) reverses the orientation. As (j) has one pole, a point w of large 
modulus has one preimage on T and the map is reversing orientation at this 
preimage. We conclude that the degree of cj) equals —1; see |26l §5] for the 
dehnition of the degree. 

Suppose hrst that all zeros of 0 are in D'^ U D~ . Equivalently, g has 
no neutral hxed points. Denote by and N~ the numbers of zeros of 0 
in and D~ respectively. Then N~^ — N~ equals the degree of 0 so that 
N~^ — N~ = — 1 and thus 

N- = N+ + 1 . ( 10 ) 

For the number N = + N~ of hxed points of in T we thus hnd that 


N = 2N^ + 1. 


( 11 ) 


Since = 1 — \dg\‘^, the zeros of cj) in are attracting hxed points of 
g while the zeros of 0 in D~ are repelling hxed points of g. Thus Lemma 
yields that 

< 2 . ( 12 ) 

It follows from ( [II| and (12) that N < 5. 

On the other hand, since 0 is odd and A-periodic it easily follows that 
that the half-periods oJi/ 2 , LZ 2/2 and lzzI 2 = (wi -f uj 2)/2 are zeros of 0. 
Equivalently, they are hxed points of g. Thus we have N > 3. Altogether, 
since N is odd by ( [II| ), it follows that N = 3 or N = 5. 

It remains to determine the criterion distinguishing the cases. Suppose 
hrst that all three half-periods are in D~-, that is, they are repelling hxed 
points of g. Then N~ > 3. This yields that N~ = 3 and A^"*" = 2 so that 
A^ = 5 by (10) and (12). Suppose now that one half-period, say 00 j/2, is not 
in D~ and thus in D^. Thus ujj/2 is an attracting hxed point of g and hence 
attracts a critical orbit by Eaton’s theorem. However, since g is odd we see 
that Uj/2 in fact attracts both critical orbits. Thus, by Eaton’s theorem, 
there are no other attracting hxed points. Thus N~^ = 1 and hence N = 3 
in this case. 

We see that A^ = 3 if and only if there exists j G {1, 2, 3} such that oJj/2 
is an attracting hxed point of g; that is, \dg{ujj/2)\ < 1. Using ([^ this takes 
the form 

p • n 

< 1. (13) 


mm 

i<i<3 
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Now 0 yields 

a rji 

b uji huji 




+ 


r^ilcjip Imr 

CJiTT 


iOl 


rjiOJi Imr ^ 


and 


ej ej IcJi p Im r uji ejOjf Im r 
b 71 Ui TT 


Substituting the last two equations in (13) yields 


ejuil + T]iUJi 
71 


Imr — 1 


< 1 . 


This is equivalent to strict inequality in ([^. This completes the proof of the 
theorem in the case that all zeros of 0 are in D'^ U D~. 

To deal with the case where this condition is not satished, we note that in 
the above arguments we may replace (f){z) by Tp^z) = (f){z) —w for any w E C. 
Noting that 0 and 0 have the same Jacobian determinant we conclude that 
whenever all tc-points of 0 are in U D~, then 0 has either three or hve 

tc-points in T. Moreover, there are three tc-points if and only if one them is 
in D~^. 

We will use the following result [3 Proposition 3]. 


Lemma 5. Let D <Z C be a domain and let f: D ^ C be a harmonic map. 
Suppose that there exists m E N such that every w E C has at most m 
preimages. Then the set of points which have m preimages is open. 


Suppose now that 0 has a zero in T\{D^VJD~). Then there are arbitrarily 
small w such that 0 has a tc-point in . Since for such w the function 0 
has three w-points, it follows from Lemma that 0 has three zeros in T. 
These zeros are the half-periods and we see that we have equality in (|^. 
This completes the proof of the theorem. 

Remark 1. The quantities CjOj'l -t- rfiOJi occurring in (|^ have the following 
representations via theta functions, see [32l P- 44]: 

2 ^2(0) 2 ^o(O) 2 ^3(0) 


The series for theta functions converge fast and provide a convenient way to 
compute Figure [Tj 
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Let 


fAr) 


ZTTZ 

ejcuj + r]iUi 


so that (|^ takes the form ImFj(r) > 0. It seems that Fj maps the compo¬ 
nents of the set of all r in the upper half-plane where Im Fj (r) > 0 univalently 
onto the upper half-plane, but we have not been able to prove this. 


Remark 2. Lin and Wang Theorems 1.6 and 1.7] pay special attention 
to the case that r = 1/2+ ib, in which case Green’s function has hve critical 
points if and only if b is outside a certain interval [bo, &i] where bo ~ 0.35 and 
6i 0.71. 

Here we note that the constants bo and bi are related to the so-called one- 
ninth constant A occurring in approximation theory; see [HI Section 4.5] for 
the definition and properties of this constant. As noted there, this constant 
was already computed by Halphen m p. 287] to six digits. It turns out 
that A = The numerical values are A = 0.10765391..., 

bo = 0.35472989 ... and = 0.70476158 .... 

To prove the above relation between these constants we note hrst that 
Ciuf -1- rjiUi = 0 if and only if '^2(0) = 0 and thus if 


^(2fc + l)2/i^(^+b = 0 where h = e*"" 

k=0 


For T = 1/2 + ib with 6 G M we have h = and hence 


j^k{k+l) 


exp ( (in — 2TTb) 


k{k + l) 




For T = 1/2 + ib the condition that eicol + rjiUJi = 0 is thus equivalent to 

00 

^(2fc + l)2(-a;)*’(*’+i)/2 ^ Q_ 
k=0 

The smallest positive solution of the last equation is the one-ninth constant A; 
see [m Section 4.5]. Since x = we deduce that A = 

The modular group leaves the sets where Green’s function has hve or 
three critical points invariant. The transformation T{z) = {z — l)/{2z — 1) 
satishes T(l/2 -|- ib) = 1/2 -|- i/[Ab). This implies that bi = l/(46o) and 
A = 
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Figure 2: Julia sets corresponding to r = i and r = 


Remark 3. Figure shows the Julia sets of the functions g corresponding to 
the square and the hexagonal lattice. The standard fundamental domains 
are marked by thick lines and the half-periods are marked by circles. For the 
square lattice, the half-periods are the only hxed points of with (a;i-l-a; 2)/2 
attracting and the two other ones repelling. For the hexagonal lattice, all 
half-periods are repelling and there are the attracting hxed points (a;i-t-a;2)/3 
and 2 {ui + 002 )/“i marked by crosses. 

Figurej^shows the Julia set of the function corresponding to r = 1/2+ihi. 
The half-period Ci;i/2 is a parabolic hxed point with two petals. The two other 
half-periods are repelling hxed points. 

Comment. Besides the application in [22[ [20l EU [19] , anti-holomorphic dy¬ 
namics were studied for its own sake in [II11301I29] and elsewhere. In these 
papers, the iteration of the anti-holomorphic map + c was investigated. 
The bifurcation diagram in the c-plane for this map is called the Mandelbar 
set or Tricorn. It was noticed that the neutral cycle appears on a whole arc 
of the boundary of the hyperbolic component. 

In the case we considered, a striking new phenomenon occurs: the neutral 
cycle occurs everywhere on the boundary of the hyperbolic component. 

To be more precise, let H be the upper half-plane, so r E H. Let X be the 
maximal open set of r where G has three critical points, and Y = H\X. For 
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Figure 3: Julia set corresponding to r = 1/2 + ibi. 


r in X, the critical orbits of g = gr tend to an attracting fixed point, while 
for r G F they tend to two attracting fixed points. The open sets X and 
Y and their common boundary which is a piecewise analytic curve exhaust 
the whole parameter space! So a generic gr in our family has hyperbolic 
dynamics in a trivial way. 

Such a situation occurs in holomorphic dynamics only in trivial cases [251 
Theorem 2.2], and was not encountered so far in anti-holomorphic dynamics. 

The point in a hyperbolic component at which the multiplier is zero is 
called the center of this hyperbolic component. Computation shows that 
the centers of X and Y in the standard fundamental domain are t = i and 
T = corresponding to the square and hexagonal lattices respectively. 

We thank C.-S. Lin for his useful comments on this paper. 
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